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Abstract. We report on the recent progress in theoretical and numerical studies of 
entanglement entropy in lattice gauge theories. It is shown that the concept of quantum 
entanglement between gauge fields in two complementary regions of space can only be 
introduced if the Hilbert space of physical states is extended in a certain way. In 
the extended Hilbert space, the entanglement entropy can be partially interpreted as 
the classical Shannon entropy of the flux of the gauge fields through the boundary 
between the two regions. Such an extension leads to a reduction procedure which 
can be easily implemented in lattice simulations by constructing lattices with special 
topology. This enables us to measure the entanglement entropy in lattice Monte-Carlo 
simulations. On the simplest example of Z2 lattice gauge theory in (2 -I- 1) dimensions 
we demonstrate the relation between entanglement entropy and the classical entropy 
of the field flux. For SU (2) lattice gauge theory in four dimensions, we flnd a signature 
of non-analytic dependence of the entanglement entropy on the size of the region. We 
also comment on the holographic interpretation of the entanglement entropy. 
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The methods of quantum information theory has led recently to many important 
advances in our understanding of quantum field theories in continuous space-times and 
on the lattices [I]. In particular, quantum entanglement of states of systems with many 
degrees of freedom turned out to be a very useful model-independent characteristic of the 
structure of the ground state of quantum fields. A commonly used measure of quantum 
entanglement of the ground state of quantum fields in [D — l)+l-dimensional space-time 
is the entropy of entanglement 5* [A] between some {D — l)-dimensional region A and 
its {D — l)-dimensional complement B, which characterizes the amount of information 
shared between A and B Entanglement entropy is defined as the usual von Neumann 
entropy for the reduced density matrix pA associated with the region A: 

S [A] = -ti A {pA^n pa) (1) 

The reduced density matrix is obtained from the density matrix of the ground state 
of the theory, pab = |0)(0| , by tracing over all degrees of freedom which are localized 
outside of A, i.e. within B [1]: 

PA = tTB PAB = tlB \0){0\ (2) 

This density matrix describes the state of quantum fields as seen by an observer who 
can only perform measurements within A. 

The entropy of entanglement of confining gauge theories has recently become a 
subject of extensive studies in the framework of AdS/CFT correspondence, where a 
simple geometric expression for the entanglement entropy was conjectured [21 [3l H]. 
One of the most interesting predictions of [21 [3l H] is that for confining gauge theories 
the entanglement entropy ([T]) should be non-analytic in the size of the region A. In 
the limit Nc —>■ oo the derivative of the entropy over the size of A changes from 
being proportional to to a quantity of order A*"^. This property can be intuitively 
understood without any reference to AdS/CFT. Indeed, at small distances an observer 
in A can only see quarks and gluons, whose number and hence the entropy scales as N^, 
where Nc is the number of colours. At large distances the effective degrees of freedom 
are quarks and hadrons, whose number and the entropy are of order N^. Thus at some 
characteristic size of A, which should be determined by a typical hadronic scale, the 
entropy should change from to Nj^. The nontrivial prediction from AdS/CFT is that 
this change is stepwise, and in some sense colourless and colourfull degrees of freedom 
never coexist at one energy scale. A similar non-analuticity has also been predicted 
using the approximate Migdal-Kadanoff decimations for SU (2) lattice gauge theory [5] 
and observed in lattice simulations [6]. Thus the behavior of entanglement entropy of 
confining gauge theories can be an interesting new test of Maldacena duality between 
string theories on {D + l)-dimensional curved spaces and Z)-dimensional gauge theories 
which live on their boundary. On the other hand, an attempt to define rigorously 
what is a quantum entanglement between gauge fields in the regions A and B also 
leads to the result which closely resembles the holographic principle as formulated by 
t'Hooft [7]. Namely, it turns out that due to a specific structure of the Hilbert space of 
gauge theories, the entanglement entropy includes the classical Shannon entropy of the 
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probability distribution of the flux of the curvature tensor of the gauge field through 
the boundary dA of the region A [8] . This entropy can be interpreted as the entropy of 
the endpoints of electric strings on dA. Thus the boundary between the regions A and 
B becomes a sort of D-brane for electric strings. 

In order to define the partial trace tr^ in ([2]), one should decompose the Hilbert 
space of the quantum field theory into a direct product Ha T^b of Hilbert spaces 
Ti-A and Tis of states of quantum fields inside the region A and B. One of our main 
assertions is that such a decomposition is impossible for a Hilbert space of physical 
states, i.e. a Hilbert space of states which satisfy the Gauss law. The reason is that 
in contrast to scalar field theories, in pure gauge theories elementary excitations are 
not associated with points in space, but rather with closed loops, which are the lines 
of electric flux or electric strings [9]. Obviously, closed loops cannot be classified as 
belonging either to A or B, and there are the states which cannot be decomposed into 
direct products of states localized completely within A or B. It should be stressed that 
this conclusion follows from the general structure of the Hilbert spaces of gauge theories 
rather than from the dynamical properties of a particular theory, such as confinement 
or deconfinement. 

For the sake of brevity let us consider the simplest case of Z2 lattice gauge theory in 
(2 + 1) space-time dimensions. All the considerations below can be generalized to other 
gauge theories with minor modifications. The Hilbert space is a space of all functions 
ijj [zi] of Z2-valued link variables zi on a two-dimensional lattice. A convenient basis 
in this space can be parameterized by a set of integer variables m/ = 0, 1 defined on 
the links of the lattice: ip [^i] ''^i] ~ 11 -^P- This basis is orthonormal with respect 

to the scalar product {ipi\4'2) = X] V^i [^i]'^2 [zi]- The Hilbert space of physical states 

7^0 is obtained by imposing the Gauss law as a first-class constraint, so that the wave 

functions of physical states are invariant under gauge transformations: ip [zp] = [zi] . 

The basis states in Ho can be labelled by such sets {mi}, for which the links with m/ = 1 

form closed non-self-intersecting loops, and can be obtained by acting with all possible 

Wilson loop operators [C] = Yl the trivial strong-coupling ground state with 

i&c 

ip [zi\ = const. The operator W [C] thus creates an electric string on the loop C. 

Let us denote the Hilbert spaces of all functions on links within A or i? as Ha 
and Hb- We would like to decompose the Hilbert space of physical states Ho as 
Hq = Ha ® Hb, so that Ha and Hb are the subspaces of Ha and Hb- Consider 

some state with the wave function of the form ip [zi\ ~ 11 -^^ where the C is some closed 

c 

loop which belongs both to A and B. This state is a direct product of two states in Ha 
and Hb with wave functions ipA [zi] ~ 11 ^^^1 ipB [zi\ ~ 11 where Ca and Cb are the 

Ca Cb 

parts of the loop C which belong to A and B. Since the basis states corresponding to 
different loops are orthogonal, this decomposition is unique, and the states ipA and ipB 
should necessarily belong to Ha and Hb- The spaces Hq and hence Ha and Hb should 
also contain the trivial strong-coupling ground state with ipo [zi] = const = ipAoi^Bo, 
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i'A,Bo [-2/] = const. Therefore if Hq is indeed a direct product of Ha and TCb, it should 
also contain a direct product of, say, ipAo with ipB- The wave function for such a state 
is proportional to Yl^h but Cb is in general not closed, therefore the Gauss law is 



violated at the boundary between A and B. We thus arrive at the contradiction, which 
completes the proof. 

Let us now try to extend the Hilbert space of physical states in some minimal way, 
so that a direct product structure can be introduced. The Hilbert space Ho cannot be 
decomposed into a direct structure, because it contains the closed strings which belong 
to both regions A and B. We can get rid of such strings, if we allow the strings to open on 
the boundary between A and B. In this case any closed string which crosses dA can be 
represented as a direct product of two open strings, each of which lies completely either 
within A or B. We thus arrive at the minimal extension Ho of Ho, which is the Hilbert 
space of states which violate the Gauss law only at dA. Minimal extension here means 
that the Gauss law does not hold in a minimal number of points, and Ho is a subspace 
of H of minimal dimensionality which contains Ho and which can be represented as a 
direct product. To see that this is indeed so, suppose that the Gauss law still holds in 
some points on dA. In this case one can repeat all the arguments above for the loop C 
which goes through these points and arrive at the same contradiction. 

We thus conclude that in order to define the spaces of states of gauge fields which 
are localized either in A or i? one should extend the Hilbert space of physical states and 
include also the states of electric strings which can open on the boundary dA between 
A and B. This boundary can be therefore considered as a sort of D-brane for electric 
strings, with the positions of their endpoints being the additional degrees of freedom 
which emerge as a result of such extension. It is reasonable to conjecture that these 
extra degrees of freedom should contribute somehow to the entropy of entanglement 
between A and B. In the extended Hilbert space one can naturally define the partial 
trace over the fields in i? as a sum over all links which belong to B: 



The same expression is valid for lattice gauge theory with any gauge group, with 
the summation over Z2-valued variables zi being replaced by integration over the link 
variables which belong to the corresponding gauge group. 

Having the reduction procedure ([3]) at hand, one can use the standard replica trick, 
which is commonly used to calculate the entanglement entropy of scalar quantum field 
theories Namely, we express the entanglement entropy in terms of the free energies 
F [A, s, T] = — \nZ [A, s, T] of the theory on a space with topology C^'^-' ® T^~^, where 
C*-'*-' is the s-sheeted Riemann surface and T^~^ is the D — 2 dimensional torus: 




(3) 




(4) 



The branching points of C^^^ should cover the boundary dA of A. 
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Figure 1. Pair production for a confining gauge theory in the vicinity of a black hole. 



In order to show that the endpoints of electric strings on the boundary of A indeed 
become physically relevant degrees of freedom, let us first consider a confining gauge 
theory - say, QCD - in the vicinity of a black hole. Since the region inside the horizon 
is not accessible to an external observer, in order to find the density matrix of the fields 
outside of the black hole one should trace over all states of fields inside the horizon, which 
is one possible way to relate the Bekenstein-Hawking entropy with the entanglement 
entropy. One possible explanation for the radiation from black holes, which is usually 
ascribed to V. Gribov, is the pair production, in which a photon emits a virtual electron- 
positron pair which becomes on-shell in the external gravitational field. One of these 
particles disappears behind the horizon and another is observed as radiated by a black 
hole. Consider now a similar process, in which a gluon emits a pair of a quark and 
an antiquark. Say, the quark disappears behind the horizon and an antiquark goes to 
an external observer. If the Hilbert space of gauge theory in the vicinity of a black 
hole were constructed in a usual way, we would have to conclude that even after the 
disappearance of the quark the particles are still connected with a confining string. 
This, however, would violate the "no-hair" theorem. There are only two ways to save 
the "no-hair" theorem: either to assume that a black hole is a source of free quarks, or 
to assume that a confining string can be torn on the horizon (see figured]). The former 
assumption is obviously prohibited by astrophysical observations, while the latter leads 
exactly to the extended Hilbert space described above. The horizon of a black hole thus 
should be a sort of D-brane for electric strings of confining gauge theories. 

Consider now a simple trial ground state wave function of (2 -|- l)-dimensional 
lattice gauge theory. This trial ground state is a superposition of all possible 

configurations of closed electric strings with the weight exp X]"^'^' where Yl^i 

is the total length of strings: 



where 6mi = means that we sum only over closed electric strings and C is the 
normalization constant. The wave function would be a product of functions of 
individual link variables if the constraint 6mi = was omitted. Thus in some sense 
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the trial ground state is only "minimally" entangled because the electric strings 
are closed, and there is no entanglement between closed string states in the bulk of the 
regions A and B. Hopefully, with such a wavefunction we can calculate the entanglement 
solely due to the strings which cross dA. 

After the application of the replica trick and some algebraic manipulations we arrive 
a the following simple and universal answer for the entanglement entropy of the trial 
ground state ([5]): 

^[^] = ~ ^ p[{Xi,...,Xm}]\np[{xi,...,Xm}] (6) 

{xi,...,Xm} 

where p [{xi, . . . , Xm}] is the probability that the electric string crosses dA in the m 
points xi, . . . ,Xm- Thus for the trial ground state ([5]) the entanglement entropy is the 
classical Shannon entropy of the string endpoints on dA. This result can be interpreted 
as follows: an observer within A can not learn whether the electric string are continuous 
over dA or not, and should assume the latter. Although in the ground state all electric 
strings are, of course, continuous, the observer who can only perform measurements 
within A can never know this, and for him the positions of string endpoints outside of 
A are an additional source of uncertainty. 

In order to check the applicability of the expression ([6]) to the true ground state 
of (2 + l)-dimensional lattice gauge theory, we have measured both the entanglement 
entropy and the entropy of string endpoints (l6l) in Monte-Carlo simulations. In order 
to measure the entanglement entropy, we have used lattices with topology C*^'*-' T'^"^ 
and approximated the derivative over s at s = 1 in (jlj) by a finite difference between 
s = 2 and s = 1. The region A was a square of size I x I, and only the differences of 
entropies S [A'] — S [A] = S {I + 1) — S (/) were measured. 

The entropy of the endpoints of electric strings was calculated as the Shannon 
entropy of the classical probability distribution p[{xi, . . . ,Xm}] of the intersection 
points Xi, . . . ,Xm between electric strings and the boundary of A under the simplifying 
assumption that at fixed m all endpoints are uniformly distributed over dA. The 
configurations of electric strings were extracted from lattice simulations of the theory 
which is related to Z2 lattice gauge theory by a Kramers- Wannier duality - namely, the 
Ising model in (2 + 1) dimensions. 

Lattice derivatives S' (Z + 1) — 5* (/) of the entanglement entropy and the entropy of 
the distribution of intersection points are plotted on figure El For the left plot on figure 
|2] the coupling constants Pg = 0.788 of Z2 lattice gauge theory is close to the critical 
point at which percolating electric strings emerge. At this value of coupling constant 
the differences S" (/ + 1) — S* (/) of both entropies indeed agree within the error range and 
practically do not depend on I, which means that the entropy is proportional to the area 
of the boundary with a good precision. We thus recover the familiar "area law" for the 
entanglement entropy. The data plotted on the right plot on figure[2]was obtained for the 
fixed size of the region A I = S and for different values of the coupling constants Pg and 
j3s. Again, only the differences S {I + 1) — S (/) were compared. The dependence on Pg is 
again similar for both entropies. We conclude that the entanglement entropy is indeed 
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Figure 2. Lattice derivatives of the entanglement entropy and the entropy of 
intersection points over the size of the region for 16^ lattice. On the left: as a function 
of I aX Pg — 0.788, on the right: as a function of (3g at I = 8. 



saturated by the uncertainty in the positions of the endpoints of open electric strings 
in the extended Hilbert space at the outer side of dA. In view of the actively discussed 
"holographic descriptions" of field theories, which relate the properties of string theories 
in the bulk of AdS space or modifications thereof with the properties of the field theories 
on its boundary, it seems rather tempting to associate the entanglement entropy of the 
region A with the classical entropy of some statistical theory on its boundary. The 
endpoints of electric strings could then be interpreted as elementary degrees of freedom 
in such a theory. 

Finally, let us present the results of the measurements of the entanglement entropy 
in four-dimensional SU (2) lattice gauge theory. As explained in the introductory 
paragraph, one can expect that for confining non-Abelian lattice gauge theories the 
entropy should change rapidly at some scale related to a typical hadronic scale. In order 
to measure the entanglement entropy in this case, we have used the same approximations 
as for the (2 + 1) Z2 lattice gauge theories considered above. 

^ ' ' a = 0.iofm > ' • 

a = 0.11 fm i -H 

a = 0.12fm > T - 

• a = 0.14fm ♦ 



T 

♦ 



.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 
I, fm 

Figure 3. The dependence of the derivative of the entanglement entropy ^'qP on 
I. Solid line is the fit of the data by the function Cl~^. 

The derivative 4 5* (Z) estimated from these measurements is plotted on figure 
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[31 It can be seen that this derivative grows rapidly at small distances. For comparison 
with the asymptotic behavior ^ 5 (/) ~ at / ^ 0, we have fitted these results by 
the function Cl~^ (solid line on figure [3]). For the data points with the smallest / the 
finite differences {F (/ + 1, 2, T) — F [/, 2, T]) were found at fixed l/a at different a, 
so that the finite differences (/ + a/2)~^ ~ ~ '^/2)~^ still behave as ~ a~^. At 
larger / ^ S"/ (Z) seem to approach a kind of plateau for the values of I between 0.3 fm 
and 0.5 fm. Here the values of ^ Sf (/) obtained for different values of lattice spacing 
differ rather significantly, which indicates that for our lattice parameters finite-volume 
and finite-spacing effects may still be rather strong. Nevertheless, at least qualitatively 
all data points for different values of a display the same behavior. At / ~ 0.5 fm the 
derivative ^ S (l) rapidly goes to zero and remains equal to zero for larger /. We 
thus find an indication of a nonanalytic behavior of entanglement entropy, in accordance 
with the predictions of [21 El IH IS] • 

The entanglement entropy is thus a new and a universal parameter which can 
be used to characterize the confining and deconfining phases of gauge theories. 
Interestingly, it should distinguish between confinement and deconfinement even in 
the presence of dynamical charges and can be therefore considered as a long-sought- 
for "order parameter" for confinement /deconfinement phase transition in QCD. On the 
other hand, the entanglement entropy receives an essential contribution from the degrees 
of freedom on the boundary of A. We can interpret this fact in terms of "holography" - 
namely, that the dynamics of gauge theories in D dimensions can be encoded on D — 2- 
dimensional manifold dA. This is at least an interesting point for further research. 
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